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1.  Introduction 


Tropical  cyclone  forecasters  ha  vs  in  the  past  relied  primarily 
on  analoq  and  statistical-climatol ogical  aids  for  objective 
guidance.  acre  recently#  dynamical  aodals  have  provided  routine 
guidance  for  dost  northern  hemisphere  tropical  cyclones. 

Development  of  these  dynamical  models  has  progressed  such  that 

dynamical  guidance  appears  to  be  comparable  or  superior  to  the 
official  foracasts,  especially  beyond  24  h  (Blsberry, 1979) . 

Southern  hemisphere  tropical  cyclone  forecast  aids  generally 
have  been  patterned  after  those  used  in  the  northern  hemisphere. 
Such  techniques  are  hampered  by  a  deficiency  of  atmospheric  data 
near  the  tropical  cyclone  and  in  the  surrounding  region.  Aircraft 
reconnaissance  has  only  recently  been  attempted  on  southern 

hemisphere  storms,  sc  that  the  locations  of  the  storm  centers  are 
less  accurate  than  for  Atlantic  and  western  north  Pacific 

tropical  cyclones. 

Among  the  current  southern  hemisphere  forecast  aids  is  an 
analog  method  used  by  the  Joint  Typhoon  Earning  Canter  called 
TYAH78  (Ocean  Data  Systems,  Inc.,  1978),  which  is  an  updated 
version  of  the  earlier  analog  schema  (Jarrell  and  Wagoner,  1973). 
Brand  and  Blelloch  (1976)  developed  a  simulated  analog  technique 
called  SWPAC  which  is  similar  to  a  climatology  and  persistence 
model  which  is  used  by  the  National  Hurricane  Cent--  for  Atlantic 
storms.  The  Australian  Bureau  of  Meteorology  has  developed  an 
analog  scheme  called  CYCLOGUE  as  part  of  its  operational  guidance 
(Annette,  1978).  Chong  et.  al.  (1980)  davised  a  regression 
scheme  which  provided  forecasts  with  smaller  errors  tnan  t in¬ 
official  Australian  Bureau  predictions.  There  are  currently  no 
operational  dynamical  model  forecasts  mad?  for  sjutnerr. 
hemisphere  tropical  cyclones. 

The  Joint  Typhoon  Warning  Center  (JTWC)  at  Guam  was  recently 
tasked  to  expand  its  area  of  forecast  responsibility  to  include 
tropical  cyclones  in  the  southern  hemisphere.  JTWC  currently  uses 


the  analog  and  statist ica  1/climatologica  1  alas  including  TYAN78 
as  quidance  in  preparing  these  forecasts.  The  success  of 
dynaaical  models  in  predicting  northern  hemisphere  tropical 
cyclone  tracks  (Harrison,  1981;  Harrison  and  Fiorino,  1982) 
indicates  a  potential  for  improving  forecasts  in  the  southern 
hemisphere  as  well. 

Durinq  the  1982-1983  storm  season,  JTWC  will  begin  operational 
testing  of  the  Navy  Nested  Tropical  Cyclone  Model  (NTCM)  on 
southern  hemisphere  storms.  This  paper  presents  an  evaluation  of 
the  NTCM  performance  for  selected  storms  in  the  Australian  reqion 
durinq  the  period  of  1975-1980.  These  results  should  provide  a 
basis  for  estimating  the  performance  of  the  model  when  it  becomes 
operational  in  the  southern  hemisphere  tropical  cyclone  basins. 
Another  objective  of  this  paper  is  to  determine  whether  the 
operationally-analyzed  wind  fields  in  the  data-sparse  southern 
hemisphere  are  adequate  tc  support  a  dynaaical  tropical  cyclone 
model.  One  might  expect  the  southern  hemisphere  storms  to  be  more 
difficult  to  predict  because  of  the  early  recurvature  and  rapid 
poleward  motion  of  many  of  these  storms.  On  the  other  hand,  on? 
expects  that  a  dynamical  model  is  more  likely  to  be  capable  of 
predicting  this  type  of  motion  than  are  the  statistical  and 
analog  techniques.  This  potential  provides  the  motivation  for 
testing  the  dynamical  model  in  the  southern  hemisphere. 

2.  The  model 

The  NTCH  is  a  three-layer,  primitive-equation  mod  si  with  a 
movinq,  two-way  interactive  fine  mesh  grid  esbeddea  in  a  coarse 
grid  channel.  Details  of  the  NTCH  are  given  by  Harrison  (1973, 
1981).  The  model  is  currently  used  for  predicting  tracks  of  north 
Pacific  Ocean  tropical  cyclones. 

To  facilitate  comparisons  with  the  northern  hemisphere 
version,  the  model  used  in  this  3tudy  was  chosen  to  be  identical 
to  the  northern  hemisphera  version,  except  for  the  trivial 
chanqes  necessary  for  the  southern  aemisphere  application.  These 
changes  include  the  direction  cf  the  bogus  verts  spi'*  d  the 
siqn  of  the  Coriolis  parameter.  Another  modiftc*  io,.  became 
necessary  due  to  the  grid  structure  of  the  tropical  analysis  that 


is  available  to  initialize  the  eodel.  The  Fleet  Numerical 
Oceanography  Center  (FNOC)  global  band  analysis  is  performed  on  a 
Mercator  projection  from  60°N  to  only  40°S.  If  the  southern 
hemisphere  application  was  simply  a  mirrcr  image  of  that  in  the 
northern  hemisphere,  all  storms  south  of  13.6°S  would  cause  the 
relocation  of  the  coarse  mesh  grid  to  extend  bayond  40°S.  As  the 
southern  edge  of  the  coarse  mesh  grid  must  be  limited  to  40°S, 
the  initial  fine  mesh  grid  position  that  is  centered  on  the  storm 
location  is  not  constrained  to  be  initially  at  the  same  coarse 
mesh  location,  as  is  the  case  in  the  northern  hemisphere 
application.  In  particular,  the  initial  position  of  the  fine 
mesh  grid  is  allowed  to  be  considerably  closer  to  the  poleward 
boundary  to  fit  the  coarse  grid  within  the  global  band  domain. 
Fortunately,  the  vast  majority  of  the  southern  hemisphere 
tropical  cyclones  experience  a  rapid  decay  well  before  they 
approach  the  southern  boundary  of  the  domain. 

The  sensitivity  of  the  model  to  a  fixed  coarse  mesh  boundary 
at  40  degrees  latitude  was  tested  for  three  northern  hemisphere 
cases.  An  example  (Fig.  1)  illustrates  that  the  predicted  track 
with  a  hypothetical  40°N  coarse  grid  limit  begins  to  depart 
significantly  from  the  standard  model  forecast  after  only  24  h  of 
integration.  The  reason  the  effects  are  felt  so  quickly  is  that 
the  application  of  a  free-slip  wall  boundary  condition  during  the 
initialization  and  during  the  integration  distorts  the  wind  and 
geopotential  fields  some  distance  from  the  poleward  edge  of  the 
domain.  Hhen  the  wall  is  placed  closer  to  the  initial  storm 
location,  a  poleward  moving  storm  comes  into  the  modified  flow 
region  earlier.  Later  in  the  integration,  the  wall  affect  becomes 
more  noticeable  as  the  boundary  condition  on  the  coarse  mesh  grid 
tends  to  inhibit  flow  normal  to  this  boundary.  As  a  consequence, 
the  storm  must  eventually  move  parallel  (either  eastward  or 
westward)  to  the  wall.  Fortunately,  this  fixed  southern  boundary 
at  40°S  will  no  longer  be  a  problem  since  the  global  band 
analysis  is  currently  being  extended  to  cover  the  entire  globe. 
One  might  expect  that  the  predictions  based  on  these  fields  will 
be  superior  to  the  tests  reported  here,  because  the  initial 
position  of  the  poleward  boundary  can  be  placed  as  far  away  as 
desired  from  the  tropical  cyclone  center. 


3.  Model  perf oraance. 

The  NT  CM  was  tested  cn  185  cases  from  34  storas  ir.  the 
Australian  region  which  occurred  during  1975-1979.  An  attempt  was 
made  to  select  cases  on  a  random  basis,  however,  the  availability 
of  archived  FNOC  data  fields  for  the  model  initialization 
dictated  to  some  extent  which  cases  were  selected. 

A.  Storms  east  of  135®E. 

Mean  error  statistics  for  the  cases  rast  of  135  E  ar  -  giver,  in 
Table  1.  The  forecast  error  is  defined  as  the  distance  between 
the  forecast  position  and  the  corresponding  best  *rrack  position. 
The  NTCM  performance  as  measured  by  forecast  error  for  these 
cases  is  not  as  good  as  in  the  northern  hemisphere  where  the 
forecast  errors  are  on  the  order  of  210,  390  and  585  kn  at  24,  48 
and  72  h,  respectively  (Harrison  and  Ficrino,  1982).  The  NTCM 
performance  appears  to  be  comparable  or  slightly  superior  to 
Brand  and  Blelioch's  (1976)  SWPAC  analog  technique  (243,  501  and 

748  km  forecast  error  at  2  4,  48  and  72  h  respectively),  and  to 

Annette's  (1978)  CYCLOGUE  analog  (272,  575  and  713  km  forecast 

error  at  24,  48  and  72  h  respectively).  It  should  be  emphasized 
that  the  error  statistics  for  these  analog  schemes  are  based  or. 
much  smaller  sample  sizes  (SWPAC  had  41,  40  and  25  cases  at  24, 
43  and  72  h  respectively;  CYCLOGUE  had  43,  22  and  21  cases  at  2u, 
48  and  72  h,  respectively).  Because  the  samples  are  different, 
one  must  excercise  caution  in  ir.terDr  eti.ag  these  -rrror 

statistics. 

In  evaluating  the  model  performance,  it  is  useful  to  consider 
not  only  the  magnitude  of  the  forecast  error,  but  also  how  -hat 
error  compares  to  the  actual  distance  the  storm  has  moved.  For 
example,  a  forecast  with  a  72  h  error  of  209  km  may  be  considered 
good  if  the  storm  has  moved  500  km  during  that  period.  However, 
the  same  200  lea  error  might  be  regarded  as  a  poor  forecast  if  “he 
storm  has  moved  only  100  km.  The  observed  distance  (Fig.  2),  a 
measure  cf  how  far  the  storm  has  moved,  is  the  iistar.es  between 
the  best  track  position  and  the  initial  storm  position.  Nctic:- 
t’nat  this  is  always  smaller  than  the  length  of  the  storm 

tralectory,  because  of  the  curved  path  between  the  two  endpoints. 
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cases  east  of  135°7 


Porecast  error  FF  (km): 
Observed  distance  on  (km) 
Ratio  FS/CD : 

Riqht  anqle  error  (ka) : 
Speed  error  (kn) : 


Number  of  cases: 


X  -  BEST  TRACK 


O  -  FORECAST  LOCATION 


Schematic  illustration  of  observed  distance,  right  ar.g 
error  and  sp-ed  error  definitions. 


B 


f 

r- 

k. 


I 

ii 
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If  the  ratio  of  the  naan  forecast  error  (?2)  to  rha  mean  observed 
distance  (OD)  is  lass  than  1,  soma  skill  in  forecasting  nay  be 
said  to  exist.  In  '■he  northern  hemisphere,  the  HTCM  exhibits 
?E/OD  values  of  .58,  .55  and  .60  at  2  4,  48  and  72  h  which  can  be 
regarded  as  a  40-45J  skill  in  forecasting.  For  these  southern 
hemisphere  cases,  the  MTCH  forecasts  have  only  a  17-20%  skill. 

The  riqht  angle  error  (Fig.  2)  is  defined  as  the  normal 
distance  froa  the  forecast  position  to  the  line  connecting  the 
initial  and  best  track  positions.  Thus,  it  is  a  measure  of  how 
well  the  model  predicts  the  direction  of  the  s*ora  action.  Th~ 
distance  alonq  this  line  froa  the  best  track  to  the  intersection 
with  the  right  angle  error  line  is  defined  as  the  sp^ed  error 
(Fig.  2)  ,  because  it  measures  the  displacement  error  that  results 
from  the  incorrect  prediction  of  storm  translation  speed.  The 
NTCM  has  a  smaller  riqht  a tqle  error  than  the  speed  arror  at  24  h 
(Table  1)  ,  but  has  a  smaller  speed  error  than  right  angle  error 
at  longer  forecast  times. 

Meal  (1977)  defined  five  categories  of  southern  hemisphere 
storm  motion:  southeastward,  southwestwari,  r=curvatura  to  th* 
southeast,  recurvature  to  the  southwest,  and  looping  motion. 
Table  2  summarizes  the  NTCB  performance  according  to  these  storm 
motion  categories.  Each  case  i3  placed  into  a  particular  category 
according  to  the  subsequent  72  h  best  track,  rather  than  its 
complete  history.  Thus,  a  forecast  case  froa  ‘■he  early  stages  of 
an  eastward  recurving  storm  which  later  fits  into  category  3  may 
first  have  been  classified  in  category  2.  The  STCB  performs  bast, 
in  terms  of  the  FE/OD  ratio,  when  the  storm  is  moving  to  the 
southwest  (category  2).  One  such  forecast  is  presented  in  Fig.  J. 
The  model  was  able  to  predict  both  kinds  of  recurving  storms 
(categories  3  and  4)  moderately  well.  The  worst  model  p-rformar.ee 
was  in  cases  of  looping  storms  (cat eg cry  5).  In  these  cases,  the 
poor  FE/OD  ratio  can  be  partially  attributed  to  the  small 
observed  distances  of  looping  storms  rather  than  to  larg- 
forecast  errors.  Storms  moving  to  *hs  scutneast  (category  1) 
were  also  handled  poorly. 

The  right  angle  and  speed  error  biases  reveal  ooma  of  the  STCa 
systematic  error  characteristics.  The  forecasts  generally  exhibit 


16 


a  poleward  bias  in  direct  ion  and  slow  translation  spsc-ds.  This 
slow  tendency  is  «or?  severe  than  for  the  northern  hemispnere 
version,  which  also  predicts  stern  direction  very  w^ii 
( Harris  on  ,1 90 1 ;  Picriro  and  Harrison,  1982)  . 

The  difference  in  performance  between  the  two  hemispheres  may 
be  due  to  the  paucity  of  atmospheric  data  available  for  the  PNOC 
southern  heaisphere  analyses.  If  data  are  sparse,  -he  resulting 
analysis  fields  are  based  aostly  on  climatology.  The  mean  summer 
southern  heaisphere  upper  atmosphere  (for  example,  Pala-n  and 
Newtcn,  1969)  is  characterized  by  a  subtropical  ridge  extending 
froa  160*E  to  17Q*H  which  produces  a  southward  or  sourhwestwar i 
environmental  flow  in  the  eastern  Australia  region.  The 
persistent  presence  of  this  cliaatol ogical  feature  in  the  model 
data  could  cause  the  poleward  track  direction  bias.  Such  a 
steering  flow  would  also  indicate  why  the  model  does  well  with 
cateqory  2  storas  but  not  with  category  1  storms. 

B.  Storas  west  of  135°E 

The  mean  error  statistics  for  the  cases  west  of  135  °  E  are 
given  in  Table  3.  The  forecast  errors  ara  smaller  at  24  h  and 
larger  at  48  h  and  72  h  compared  to  the  eastern  region  casis 
(Table  1).  Notice  also  that  the  forecast  error /observed  iistar.es 
ratio  shows  coaparabla  skill  at  24  h  but  has  values  exceeding  1.0 
at  48  a  and  72  h.  The  NTC3  did  not  perform  as  well  as  the  3WPAC 
analog  (179,  389,  and  595  ka  forecast  error  at  24,  43  and  72  h 

for  92,  92  and  66  cases)  (Brand  and  3lellocb,  1976),  or  the 

CYCLOGUE  analog  at  forecast  intervals  other  than  an  24  a  (276, 

481  and  623  ka  error  at  24  ,48  and  72  h  for  74,  55  and  32  casts) 

(Annette,  1978).  As  pointed  out  by  Neal  (1  977),  73".  of  'hi 

storas  in  ths  western  Australia  region  move  to  the  southwest,  or 
recurve  to  the  east  versus  only  39  %  of  those  in  the  eastern 

reqion.  Thus,  the  storm  tracks  ia  -he  western  area  are  jior; 

uniform  in  behavior  which  is  an  advantage  for  analog  methods  in 
this  region. 

It  can  be  seen  from  the  right  angle  and  sp-el  errors  that  th* 
NTCH  model  predicts  speed  of  translation  better  than  direction. 
Compared  to  the  eastern  region  cas-»s,  the  decrees-  in  forecast 
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Tabla  3.  Jlean  error 
cases  west  of  135®2. 


Forecast  error  F3  < Hia )  ; 
Observed  distance  OD  (km) 
Ratio  F2/0D: 

Riqht  angle  error  (km)  : 
Speed  error  (km): 


Number  of  cas  ?s 


error  at  24  h  is  associated  with  an  improved  speed  error,  and  t he 
poorer  performance  at  later  forecast  tiaes  is  due  to  at.  increase! 
riqht  angle  error. 

Almost  halt  of  the  cases  in  the  western  region  sample  we  r- 
storas  which  skirted  alonq  the  northwest  coast  of  Australia.  One 
such  case  is  presented  in  Fig.  4.  The  NTCH  has  a  ?cry  strong  bias 
to  the  left  in  these  cases,  which  accounts  for  much  of  th* 
forecast  error  in  the  sasple.  The  poor  performance  of  the  NTCH  ir. 
these  cases  is  somewhat  surprising,  since  the  NTCS  handled  this 
track  type  (Table  2,  category  2)  very  well  for  the  storms  east  of 
Australia.  If  the  steering  action  in  the  model  data  is  correct 
in  these  cases,  then  onshore  action  of  these  cases  should  result. 
One  would  expect  rapid  decay  of  the  storms  in  such  events.  In  the 
actual  cases,  the  storas  do  not  move  onshore,  apparently  because 
that  portion  of  the  circulation  which  remains  ov^r  the  warm  ocean 
heat  source  is  sustained.  Thus,  there  may  be  a  propagation  of  the 
circulation  alonq  the  coast  due  to  continued  redevelopment  ov:-c 
this  area.  The  NTCH  would  be  unable  to  simulate  this  process 
since  it  does  no*  include  such  differential  heating  due  tc 
land/ocean  effects. 

C.  Comparison  with  TYAN78 

The  Joint  Typhoon  Warning  Center  ir.  Guam  ha?  used  ar.  anaic 
method  named  TYFOON  since  1970.  The  scheme  has  sire®  undergone 
several  updates  and  modifications,  the  most  recent  cf  which  is 
called  7YAN78.  The  scheme  currently  provides  guidance  for  -he 
northwest  Pacific,  northeast  Pacific,  southwest  Pacific, 
southwest  Indian  and  north  Indian  Ocean  regions.  Records  of 
TYAN78  performance  in  the  southwest  Pacific  ar?  no*  avaiiacl?. 
Secent  performance  of  the  analog  in  the  northwest  Pacific  region 
indicates  average  forecast  errors  of  240,  450  and  685  *a  at  2u, 

48  and  72  h,  respectively  (Annual  Typhoon  Report,  1979;  Annual 
Tropical  Cyclone  Report,  1980-1981). 

To  evaluate  NTCH  performance  v-rsus  TYAN75  in  tht  southwest 
Pacific  region,  the  analog  technique  was  run  for  the  cases  in  th- 
NTC8  sample.  Because  TYAN78  requires  48  h  of  nr-vious  positional 
data,  some  cases  coull  not  be  run.  These  car?  is  uvs  n*  ?r  resuv.  i 
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from  the  NTCM  sample,  as  were  cases  for  which  analogs  could  no: 
be  found  on  the  record  tapes.  Thus,  the  following  comparisons  are 
for  a  homogeneous  sample. 

The  mean  error  statistics  for  the  HTCJ1  and  rYAN^a  cases  east 
of  135°S  are  presented  in  Table  4.  The  number  of  cases  available 
for  comparison  is  considerably  smaller  than  for  the  oriqinal 
sample,  especially  at  72  n.  The  NTC1  error  statistics  are  similar 
to  those  of  the  total  sample  {Table  1)  despite  the  smaller  sample 
size.  The  analog  method  performs  better  in  terms  of  the  for  -cast 
error  and  7E/0D  ratio  than  the  NTCtt  at  24  h,  about  the  sans  at  u8 
h,  and  worse  at  72  h  (Table  4).  Thus,  the  analog  method  performs 
progressively  worse  as  the  forecasts  move  farther  from  the 
initial  and  48  h  history  positions  upon  which  they  are  based. 

The  mean  STCH  and  TYAH78  error  statistics  for  a  homogeneous 
set  of  cases  west  of  135°B  are  presented  in  Table  5.  Again,  the 
number  of  cases  in  the  homogeneous  sample  is  much  small ar  than  in 
Table  3.  The  NTCM  error  statistics  are  somewhat  degraded 
relative  to  those  of  the  tctal  sample  (Table  3)  especially  at  48 
and  72  h.  Furthermore,  the  TYAN73  results  ace  much  better  than 
the  8TC3  in  terms  of  the  forecast  error  and  FZ/no  ratios  at  all 
forecast  times  (Table  5) .  Thus,  the  TYAU73  performs  as  well  on 
the  western  region  cases  as  did  the  Cl CL3G0E  and  SWPAC  analog 
techniques. 

The  southern  hemisphere  TYAS7S  performs  much  bett-c  or.  *  h 
western  reqion  cases  than  it  generally  did  in  -he  northern 
hemisphere.  This  may  be  d  urn  to  the  more  uniform  behavior  cf  me 
western  region  cases  as  mentioned  earlier,  or  else  tha-  th 
sample  is  too  small  for  valid  comparisons. 

D.  Comparison  with  TCWC  official  forecasts. 

Records  of  official  forecasts  made  by  tha  Australian  3ureuu  o 
Heteroloqy  Prcpical  Cyclone  warning  Centres  (?C»C)  are  no 
readily  available.  Annette  (1978)  included  official  and  CYCLDCJ" 
forecast  error  s*ati3*ics  for  five  storms  during  1976-1377.  iabl.- 
6  summarizes  the  NTCH  prefcrmance  on  these  storms  (minus  tropical 
cyclone  Ted,  for  which  PNOC  data  were  unavailable)  v-rsus  the 
TC»C  official  forecasts  and  CYCL03UE  forecasts. 
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Table  4.  Mean  error  statistic 3  for  a  hociogeneous  sa.d 
southern  hemisphere  'JTC'l  and  TYAM71  cases  east  of  135®  z. 
statistics  are  to  the  left  of  the  slash;  ?YA'J73  -o  the  riqli 

24h  4?h  72h 

247/2  06  452/44  3  679/747 

315  611  375 

.7  4/.  62  .74/.  7  3  .  73/.  35 

14  VI  39  310/3  3u  4  39/57  9 

173/1  29  272/256  436/403 

6  9  5’j 


Forecast  error  Ft  (km) ; 
Observed  distance  CD  (km)  : 
Ratio  FH/OD; 

Bight  angle  error  (kin)  ; 
Speed  error  (km): 


Number  f  case: 


73 


1  J  0 


Tabla  !5;  Similar  to  Table  4 


Forecast  error  FE  (’em): 
Observed  distance  OD  (k-n)  : 
Hatio  FE/CD: 
aiqht  angle  error  (ka) : 
Speed  error  (km): 


exc -nr  t  or  ensos  j  ■;s*.  o: 


2’1  h 

4*h 

-T  1  V. 
c.  .. 

213/1 97 

54  3/3  ^ 

7  6  7/62  6 

273 

5  35 

5  17 

.78/.  72 

1. 31/. 74 

.  34/.  7  7 

14, -i/1  25 

462/275 

n  ao/4o  4 

125/1 25 

195/2 35 

231/3*9 

4  2 

3  * 

o  j 

Number  of  cases 


Table  6.  Comparison  of  mean  forecast  errors  (km)  for  4  storms  from  1977 
(number  of  forecasts  in  parentheses).  Statistics  for  the 


The  MTCM  does  as  well  as  the  analog  technique  an  12  h. 
However,  CYCLOGUE  performs  progressively  better  a*  Ian -r  forecast 
times,  with  especially  noteworthy  performance  at  36  h  for  ffar-r. 
and  Leo.  It  should  b=  noticed  that  all  of  these  cis-s  are  in  tae 
proximity  of  the  coast.  Ivan  so,  the  HTCH  forecast  errors  are 
about  the  sane  as  tie  official  forecast  errors  at  12  and  24  a. 

4.  Post-Processinc  of  Model  Forecasts 

Elsberry  and  Frill  (1980)  proposed  a  statistical  -echr.ique  for 
post-processing  the  tropical  cyclone  tracks  of  the  FMOC  Tropical 
Cyclone  Modal  (TCM) .  The  technique  uses  multiple  linear 
regression  equations  to  remove  systematic  bias  in  whe  TCM  track 
forecasts.  The  predictors  used  in  the  equations  ar  -■  the  storms' 
initial  latitude  and  longitud®  plus  various  zonal  and  meridional 
components  of  -he  mode  1- predicted  storm  displacement  and 
velocity.  The  most  important  predictors  are  these  which  result 
from  comparison  cf  the  best  track  with  the  forecast  found  by 
integration  of  the  model  backward  in  time. 

Peak  and  Elsberry  (198  1)  applied  the  same  technique  tc  the 
northern  hemisphere  version  of  the  NT  CM,  but  with  backward 
extrapolation  positions  instead  of  backward  integration 
positions.  The  technique  decreased  the  NTCM  72  h  forecast  error 
by  110km.  If  systematic  bias  does  exist  in  the  southern 
hemisphere  tracks,  the  same  technique  should  improve  the 
forecasts. 

A.  East  of  135°E. 

The  error  bias  characteristics  cf  the  southern  hemisphere  :<7CM 
forecasts  east  of  135#E  are  shown  in  Table  7.  The  number  of  cas  ;-s 
is  slightly  smaller  because  those  cases  without  a  36  h  history 
are  inappropriate  for  the  postprocessing  scheme  and  wer-  removed. 
Thus,  the  forecast  error  for  these  cases  is  also  slightly 
different  from  these  of  the  totai  sample  (Table  1)  .  Th  ;•  zonal 
error  bias  (Table  7)  indicates  that  -he  forecasts  have  a  strong 
bias  to  the  east  of  the  best  track,  especially  a*  72  h.  The 
meridional  error  bias  indicates  a  slight  southwar  1  bias  at  3o  and 
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Table  7.  *aar.s  (x)  an  1  standard  lev  iatior.s  (&)  (<i)  ~>-  icrecas 
error  and  zonal  and  aeridior.al  -rror  biases  a 
henisphere  NTCM  cnsjs  aas-  of  135*F. 
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48  h  and  a  northward  bias  at  72  h.  This  is  similar  to  the  bias  of 
the  northern  hemisphere  version,  although  the  zonal  errors  are 
larqer,  as  is  the  72  h  meridional  bias. 

The  110  cases  are  randomly  divided  into  a  7  3-case  dependent 
sample  and  a  37-case  independent  sample.  It  can  be  s3er.  in  Fig.  5 
that  the  systematic  zonal  and  meridional  bias  trends  are  similar 
between  the  samples,  although  the  zonal  error  is  larger  in  the 
dependent  sample  and  the  36-43  h  meridional  errors  of  the 
independent  sample  are  positive.  These  differences  may  indicate 
the  necessity  of  a  larger  sample. 

Zonal  and  meridional  regression  equations  were  derived  for  the 
dependent  sample.  The  reduction  ia  variance  by  the  regression 
equations  ranged  from  29%  to  62%  and  averaged  473.  This  reduction 
in  variance  is  considerably  greater  than  that  of  the  regression 
equations  derived  for  the  northern  hemisphere  version  (?ea<  and 
Elsberry,  1981).  The  regression  scheme  generally  reduces  both 
the  means  and  standard  deviations  of  the  zonal  and  meridional 
error  biases  in  the  independent  sample  (Table  8).  Some  slight 
increases  in  bias  result  from  the  different  bias  characteristics 
of  the  samples. 

The  mean  STCS  forecast  errors  of  the  indsoen  lent  sample  with 
and  without  the  regression  adjustment  are  about  the  same  as  the 
comparable  errors  from  the  dependent  sample  (Table  9).  The  post¬ 
processing  decreases  the  aean  forecast  error  of  the  dependent 
sample  by  63  km  at  12  h  to  a  maximum  of  152  km  at  72  h.  The  same 
error  decrease  occurs  on  the  independent  sample  except  for  a 
smaller  122  km  decrease  at  72  h.  Thus,  the  NTCS  forecast  errors 
with  statistical  post-processing  are  considerably  less  t.ian  for 
the  TYAN78,  SWPAC  and  CYCL0G02  analogs.  Th-:  improvement  in 
forecast  error  reduces  the  F3/OD  ratio  to  ,o3,  .52  and  .69  at  24, 
48  and  72  h.  An  example  of  the  forecast  iapro  v-men4-  made  by  the 
regression  scheme  is  presented  in  Fig.  6. 

Scatterplots  of  the  unmodified  independent  sample  forecast 
errors  (Fig.  7)  reveal  that  the  regression  scheme  improves  th- 
forecasts  of  almost  75%  of  “he  cas-s  at  24,  43  and  72  h. 

Different  bias  characteristics  between  the  dependent  and 
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independent  samples  accounts  for  the  slight  degradation  of  25".  of 

the  cases. 

B.  West  of  13  5°E. 


The  error  bias  character istics  of  the  applicable  southern 
hemisphere  NTCH  forecasts  west  of  135°  2  are  shown  in  Taaie  10. 
There  is  an  eastward  bias  froa  24-43  h  and  a  vary  iarg-  westward 
bias  at  72  h.  The  meridional  bias  indicat  as  a  v-ry  i arg 
southward  bias  that  increases  with  tin?.  Cases  with  such  strong 
bias  should  be  iaproved  by  the  postprocessing;  however,  th-  small 
sample  size  in  the  western  region  nay  hinder  the  performance  of 
the  scheme  with  independent  data. 

The  52  cases  are  divided  into  a  35  case  dependent  sample  and  a 
17  case  independent  sample.  The  systematic  zonal  and  meridional 
bias  trends  of  the  samples  (Fig.  3)  are  similar  up  to  48  h.  There 
is  a  smaller  independent  sample  meridional  bias  at  60  and  72  h 
and  a  smaller  72  h  zonal  bias.  The  60  h  zonal  biases  are  of 
opposite  sign.  These  differences  are  probably  due  tc  tne  small 
sample  size. 

Zonal  and  meridional  regression  equations  wer-  derived  for  th:- 
dependsnt  sample.  The  reduction  in  variance  by  the  equations 
ranged  from  26%  to  90%  and  averaged  57%.  This  reauction  in 
variance  is  larger  than  that  for  -he  eastern  region  cases, 
possibly  due  to  the  extreme  systsmatio  bias  in  the  meridional 
direction . 

The  regression  scheme  is  generally  successful  at  reducing  th; 
means  and  standard  deviations  of  the  independent  sample  oinsrs 
(Table  11).  The  scheme  performs  poorly  whir;  there  is  a  bias 
difference  between  the  dependent  and  independent  samples,  waich 
again  indicates  the  need  for  a  larger  sample. 


The  mean  forecast  errors  of  the  independent  sample  are  ubcu* 
the  same  as  for  the  dependent  sample  from  12  to  36  1  but  grow 
considerably  smaller  from  48  to  12  h  (Table  12).  The  post¬ 
processing  decreases  the  forecast  error  in  the  dependent  sample 
by  40km  at  1 2  h  up  to  5  60km  at  72  h.  The  independent  simple 


Tabla  10.  Similar  to  Table  7  exc=o 
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Time  (h)  of  Cases 
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Table  11.  Similar  to  Table  H  exceo 


Forecast  Number 
Time(h)  of  Cases 


\'TC‘< 

Zonal  :*sr 

T  <r 


forecast  error  are  reductions  are  431cm  at  12  h  and  increasing  to 
167km  at  48  h.  The  error  reduction  then  decreases  to  13Bkm  at  72 
h.  The  forecast  errors  are  reduced  by  the  pest-processing  to 
values  smaller  than  the  TYAN78,  SWPAC  and  CYCL03U E  analogues.  The 
improvement  in  forecast  error  reduces  the  PE/DD  ratio  to  .40,  .58 
and  .73  at  24,  48  and  72  h.  An  example  of  ■‘•he  forecast 
improvement  made  by  the  regression  scheme  is  in  Fig.  9. 

Scatterplots  of  the  unmodified  independent  sample  forecast 
errors  versus  the  regression  modified  forecast  errors  (Fig.  10) 
indicate  that  75X  of  the  regression  modified  cases  have  smaller 
forecast  errors  than  the  ’unmodified  cases  at  24,  43  and  72  h. 
Unfortunately ,  four  cases  of  excellent  72  h  NTCH  forecasts  are 
degraded  by  the  application  of  the  regression  adjustment. 

5.  Summary  and  Conclusions 


The  performance  of  the  Navy  Nested  Tropical  Cyclone  model  is 
evaluated  for  southern  hemisphere  tropical  cyclones.  For  storms 
east  of  Australia,  the  NTCS  mean  forecast  error  is  generally  less 
than  the  SWPAC  and  CYCL0GU3  analog  methods.  In  a  homogeneous 
comparison  with  TYAN78,  the  NTCS  tracks  were  worse  at  early 
forecast  times  and  better  at  late  forecast  times.  *'«st  of 
Australia,  the  model  generally  performed  worse  than  the  analog 
techniques.  In  a  limited  comparison  of  12  ana  24  h  forecasts, 
the  NTCH  had  errors  comparable  to  the  TCNC  official  forecasts. 


The  NTCH  tends  to  have  a  poleward  directional  bias 
predicted  tracks.  This  bias  may  be  attributed  to  the 
observations,  which  causes  the  analysis  scarce  to  re 
climatoloqica  1  values.  The  NTCS  also  dil  not  forec.as 
tracks  well  near  the  Australia  coast,  especially  in  the 
cases,  presumably  due  to  lack  of  considarati or.  of 
effects  in  the  model.  It  appears  from  these  results 
NTCH,  despite  the  paucity  of  atmospheric  lata,  lo»s  s 
skill  in  predicting  tracks  of  sou-hern  hemisphere 
cyclones,  but  net  to  -he  extent  observed  in  tiw 
hemisphere. 
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Statistical  post-processing  of  tha  model  forecasts  resul-s  ir. 
a  significant  reduction  in  mean  forecast  errors.  In  the  region 
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east  cf  135  B  the  mean  errors  are  improved  by  as  much  a?  150km  at 
72  h.  The  *  ester  r.  region  forecast  improvement  is  ever,  greater, 
such  that  the  regression  modified  NT  CM  forecasts  are  superior  to 
the  analcq  forecasts. 

It  must  be  stressed  that  this  evaluation  is  preliminary  and 
can  only  be  used  as  an  indicator  of  the  future  operational 
performance  of  the  model.  A  number  of  factors,  including  the  use 
of  less  accurate  warning  track  fixes  and  the  use  of  a  new  global 
band  analysis  scheme  will  have  a  significant  impact  cr  the  real¬ 
time  NTCM  forecasts.  Nevertheless,  the  results  presented  here 
are  ar.  encouraging  indication  that  the  operational  implementation 
of  the  NTCH  will  lead  to  improved  southern  aemisphera  tropical 
cyclone  warnings. 
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